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NOTES ON THE THEORIES OF JUPITER AND SATURN. 



BY G. W. HILL, NAUTICAL ALMANAC OFFICE, WASHINGTON, D. C. 

On account of their large masses and the near approach to commensu- 
rability of their mean motions, Jupiter and Saturn offer the most interesting, 
as well as the most difficult, field for research, in the planetary perturbations 
of the solar system. In the following remarks, without treating the subject 
in a complete manner, which would be impossible here, I intend only to 
point out a method of procedure and give a few illustrations of its use. 

At present we shall notice ordy the introaction of the sun, Jupiter and 
Saturn. It will facilitate matters much if we employ differential equations 
in which the potential function is the same for both planets. This is ac- 
complished by an orthogonal transformation of variables. Let us suppose 
that the coordinates of the sun in space are denoted by 

X, Y and Z, 
those of Jupiter by 

X + X -{- xx', Y + y + xy' and Z + z + xz', 

and those of Saturn by 

X + x'+ XX , Y + y'+ xy and Z + z'+ xz , 

where z is a small constant to be so determined that the variables x, x', &c., 
may be orthogonal. 

M, m and m' denoting severally the masses of the sun, Jupiter and Sat- 
urn, the vis viva T of the system is represented l>y the equation 

2Tdt = MdX^+m{dX+dx-\-xdx'Y^m'{dX+dx'-txdxf 

-\- similar terms in F, y, y' and Z, z, ?' 
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-\- similar terms in Y, y, y', Z, z, z' . 

In order that tiie system of variables may be orthogonal, the coefficient 
ot dxdx' in this expression must vanish, which gives us, for the determina- 
tion of X, the qiiadratic equation 

'M , M 
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Of this the smaller root must be taken. Employing Bessel's values of the 

— = 1047.879, ~ = 
m m' 

4551.479X + 1=0, 



M M 

masses of Jupiter and Saturn, — = 1047.879, ~ = 3501.6. Hence 

mm' 



whence we get 

X = ,,,,^,^„ + ( ,„! ,^„ ) % . . . = 0.0002197088. 
4551.479 \ 4551.479/ 

For brevity we will put 

{m-\-xm')^ 



H = m -\- x^m' — 



M-Ym-^-m" 



I 1,2 {m'-\-xm)^ 

ju'= m' + x'm — ^ ' 



M-\-m-\-m'' 

When the numerical values are substituted these equations give 

H = 0.9990467623m, //' = 0.9997145123m'. 

As we do not wish to know X, Y and Z, but only the six variables x, y, 
z, x', y' and z' which assign the position of Jupiter and Saturn relatively to 
the sun, we can altogether neglect the first term in the last expression for 
T, and write 

rj,_ dx^^-\-df-\-dz^ , , dx"'-^dxi"' +dz''' 
^ 2df ""^^ 2de~' 

If we put 

^+y^+z^ = r^, x'^-j-y'^-^z'^ = r'^, xx'-\-yy'-\-zz' = rr's, 

the expression of the potential function is 
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Q Mm , Mm' .mm' 1 

\r''-\-2x7-r's-j-xV'Y' [^'^+ '■^''-rr's +xV^]'^ 1 —x [r"^ — 2rr's-|-r^]*^* 

From the last term it will be seen that the motion of two planets, whose 
coordinates are severally x, y, z and x' , y', z', relatively to each other, is ho- 
mothetic with the relative motion of Jupiter and Saturn. The differential 
equations of motion are 

<Py ^ IdJl 
dt' fx dy' 

dhj ^ IdJl 
dt^ f/dy' 

The first two radicals in i? may be expanded in series proceeding accord- 
ing to ascending powers of z; and, since this constant is so small, the cube 
and higher powers of it may be neglected. Thus 

g = J^.+ i^ _ xMm rLs- xMm' ^s + ^x^m ''-^ (l+3s^) 

+ ly?Mm' 43(H-3s') + ^ , ., .-} , , ,-,. ■ 

If for i? are substituted only the first two terms of this expression, the 
differential equations are easily integrated, and the variables x, y, z and x', 
y', z' represent the motion of two planets moving according to the laws of 
elliptic motion, whose mean motions are 
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In terms of symbols whose meaning is well known, we will put 
L = T/[ifm/ia], L' = i/lMm'ju'a'}, 

G = i/lMmfxa{l—e')l Q' = y'lMm'fx'a'{l—e")], 

H= i/lMmfxa{l-e^)']cosi, H' = ,/[JlfOTya'(l— e'')]cosi, 
and denote the mean anomalies by I and I', the distances of the perihelia 
from the nodes by g and g' and the longitudes of the nodes by h and h', 
and, moreover, put 

, mm' 1 



1 — X \r"^ — 2rr's + r^]^' 

We have then the following system of differential eejuations for determin- 
ing the elements L, Q, H, L', G', H', I, g, h, V , g', h' : — 
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in which it is understood that B is expressed in terms of the elements L, 
G, H, &c. 

As r is a function of the three elements L, G, I only, and r' of L', G', I' 
only, it follows that the six elements H, g, h, H', g' and h' enter in B only 
through s; hence we have the equations 

dB^^dBds^ dB _ dB ds 

dg ds dg' dH ds dH' 

dB dB ds dB _ dB ds 

'd^'~~d^d^" dH'~"didH" 

The expression for s being given by 

rr's — xx' -\- yy' -\- zz' , 
and V and v' denoting the true anomalies, the rectangular coordinates have 
the equivalents 

X = r[cos h cos (»+^) — cos i sin h sin ('y+5')], 
y = r[sin Acos(«-|-^)-|-cosicosAsin {v-\-g)\ 
2 =rsinz sin(»+5')) 

aj/-— ).'[cos A'cos {v -\-g') — cos i'sin A'sin (»'+.9'')]) 
y'= r'[sin /(.'cos («'+^')+cosi'cos A'sin (y -\-g')'\, 
z' = r'sin V sin {v -\-g'). 
Whence the following expression for s : 

s = cos(/i — K)co&[v-]-g)cos{v'-\-g') -f cos i cos i'cos{h — h')siv{v-i-g)sm{v'-\-g') 
+COS i sin (h — h') cos {v-\-g) sin {v'+g) 
— cos i sin (A — h') sin {v-\-g) cos {v'-\-g') 
-|- sin i sin i'sin {v-pg) sin {v'-\-g'). 
Kemembering that v and v' contain only the same elements as r and r', 
and that 

cos i = f , sin i = ^(g;-^A cos i' = f ;, sin i' = l/(^J=^), 
it will be found that 

^r i/ {G'~H^ j cos h-j-y{G"—H"} cos h'l = 0, 
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dt 



ri/{G'—H'}smh+i/{C"-H'')smh'~\ = 0, 



= 0. 

Hence we have the following integrals of the differential equations, 
!/( G^ — H^) cos h-\-i/{ O'^ — H'^) cos A = a constant, 
l/(G^—H^) smh+-i/{G'^—H'^)sm A = a constant, 

II-\-H' = a constant. 
These integrals may be employed to diminish the number of differential 
equations. Thus far the system of planes, to which x, y, z, &c. are referred, 
has been left indefinite: let us now assume that the plane of maximum 
areas, called by Laplace the invariable plane, is chosen for the plane of xy. 
In this case it is well known that the constants of the first two of the inte- 
grals, given above, become zero. Then we shall have 

l/(G^— ii^)cos A + v/((?"— J?''jcosA' = 0, 
/((?2— iJ2)sinA + i/(G'='— iJ'2)sinA' = 0, 

fi" + S' = c, 
c being an arbitrary constant. But since i and % are supposed contained 
between 0° and 180°, tlie radicals in these expressions must be taken pos- 
itively. Consequently the equations are equivalent to 

A' = A -F 180°, JT + If' = c, H—H'=: ^'—^'\ 

c 

These equations determine the values of the elements H, H' and A' in 
terms of the rest, and they may be used to eliminate them from R. Then 
it is plain, from the expression of s, given above, that A will also disappear 
from R, and we shall have 

R = function (i, G, L', G', I, g, t, g), 
and s takes the much simpler form 

s = —COS {v—v +g—g) + - '^g^ sin (v+g) sin {v -{-g ). 

As to the partial derivatives of R with respect to L, L', I, I', g, g', they 
are evidently unchanged by this elimination of the elements H, H', A, h'. 

But (-Tp) and (-r™] denoting the derivatives of i2 on the supposition of 

its containing the elements H, H', h, A', we have 

(dR\dR_dRdH_ dR dH' 
\dGI dG dHdG dH' dG' 
(dR\ _dR dR dH _ dR dH' 
XdG'l dG' dHdG' dH' dG' ' 
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But we also have 

dR _dB^ _ d{h'—h) ^ Q 



dH dH' dt 



hence 



(dR\ _ dB _ dB d{H+H') ^ dR^ 
[dGl dG dH da dG' 

(dB \ _ dB _ dR d jH+H') _ dB 
dG' dH dG' ' dG" 



/dB\ 
KdG'l 



Moreover 



dR _ dB dH dB dH _ dB d{H+H' ) _ dB 
do dH do dH do dH do dH' 

Thus the system of diiferential equations still retains its canonical form, 
and is 

dL „ dB dU _ dB dG __ dB dG ._ dB 
dt dl ' dt dt' ' dt dg' dt dg' 

dl ^ _dR d[_ ^ _dR dg ^ __dB^ dg_ ^ _dB 
dt dL' dt dL" dt dG' dt ' dG" 

After this system of eight differential equations is integrated, the value of 
h is found by a quadrature from the equation 

dh ^ dB 
dt do ' 

These integrations introduce nine arbitrary constants which, together with 
c, make ten. The reference of tlie coordinates to any arbitrary planes 
introduces three more, but one of these coalesces with the constant which 
completes the value of h. 

The time t does not explicitly enter B, hence the complete derivative of 
it with respect to t is 

dB _ d^dL_ dB_d]^, ^^ 
dt dL dt '^ dl dt 

If, in this, are substituted the values of -— -, -^, &c. . , from the equa- 

(Xt (Xv 

tions just given, we shall find that it vanishes; hence 

B = a constant 

is an integral of the system of differential equations. This integral may be 
employed to eliminate one of the elements, as L, from the equations. We 
can also take one of the elements, as I, for the independent variable in place 
oft. The system of equations, to be integrated, is then reduced to the six 
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A simpler form can be given to them. If the solution of i2= a constant 
gives 

L = function {L', G, G', I', g, g', I), 

and L is supposed to stand for the right member of this, the foregoing equa- 
tions can be written 
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When the values of L', G, G', V, g and g' in terms of I have been de- 
rived from the integrals of these, they can be substituted in the equation 

-J- = -j—, which will then give t in terms of ^, by a quadrature. By in- 

verting this we shall have I in terms of i; and by substituting this in equa- 
tions previously obtained we shall have the values of all the other elements 
in terms of <. 

It will be noticed that i? is a homogeneous function of i, L',G, G' and 
c of the degree — 2; hence we shall have 

r dB , r'dB , ridB , ^.dR . dR „„ 

^dZ+^dL'+^dG + ^dG'-^'-d^^-^^^^'''''''- 

and, a? a consequence of this, 

L^+L'4^+G^+G'-% + o^^ 2i2 = a const. 
dt dt dt dt dt 

Thus, if the rate of motion of each elementary argument I, I' &c., be mul- 
tiplied by the element which is conjugate to it, the sum of the products is 
invariable. 

The sines of half the inclinations of the orbits on the plane of maximum 
areas are 
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Thus, in the special case where the two planets move in the same plane, 
we have 

This equation may be employed to eliminate one of the elements Q or G' 
from B. In the same case, the expression for s is reduced to 

s = — cos(?; — v' -j- g — g'). 
Then, if we put 

G-G' = I, g-9'^h 

R will be a function of L, L', F, I, I', y, and we shall have, for determin- 
ing these variables, the system of diiferential equations 

dL _ dR^ dL___ ^ c^/" ^ dR 

~dt lil' "dt llT' dt ~~ df ' 

dl dR dl' dR dy dR 

~dt ~dL' Tt liL" ~dt " df' 

After these ai-e integrated, the value of ^+g'' will be got by a quadrature 
from the equation 

%+^ _ __dR 
dt do ' 

If the value of L is obtained from the solution of i2 = a constant, and 

we have 

L = function {L', F, I', y, I), 

and I is adopted as the independent variable in place of t, the solution of 
this special case is reduced to the integration of the four equations 

dL' dL dl' dL dP dL dy dL 

~dl~'dT' lU" dJ7' ~dT~'df' 'dl~ dT' 
The angle between the planes of the orbits of Jupiter and Saturn is about 
1J°. This is small enough to make the terms, which are multiplied by the 
square of the sine of half it, and which are besides of two or more dimen- 
sions with respect to disturbing forces, practically insignificant. Thus, while 
we are engaged in developing those terms of the coordinates which demand 
the highest degree of approximation relatively to disturbing forces, we shall 
assume that the planes coincide; the determination of the effect of non-co- 
incidence of these planes being reserved to the end, when it will be always 
sufficient to limit ourselves to the first power of the disturbing force. 

[To be continued.] 



